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Abstract.  Detailed  models  of  the  crust  and  upper  mantle  are  crucial  for  precise  event 
location  and  computation  of  realistic  synthetic  seismograms  for  earthquake/explosion 
identification.  However,  such  models  are  not  available  for  most  of  the  globe,  including 
eastern  China.  In  this  report,  data  from  the  1992  Yellow  Sea  earthquake  is  used  to  put 
constraints  on  the  upper  mantle  velocity  beneath  the  Yellow  Sea  and  eastern  China.  The 
average  velocity  in  the  upper  mantle  is  measured  at  8.06  km/sec.  Comparison  of  observed 
waveforms  to  synthetics  suggests  that  the  upper  mantle  P-wave  velocity  gradient  is  0.0013 
kmisec/km,  with  a  higher  gradient  in  the  uppermost  mantle  just  beneath  the  Moho. 
Availability  of  array  data  for  a  suite  of  events  could  be  used  to  place  strong  constraints  on  the 
velocity  structure. 
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OBJECTIVE:  To  determine  the  crustal  and  upper  mantle  structure  beneath  eastern  China 
PRELIMINARY  RESEARCH  RESULTS: 


Introduction 

The  details  of  the  crust  and  upper  mantle  structure  of  eastern  China  have  been 
measured  in  only  a  few  places  (Mangino  and  Ebel,  1992  and  references  cited  therein). 
Mangino  and  Etel  (1992)  used  teleseismic  receiver  function  observations  to  estimate  the 
shear  velocity  of  the  crust  at  stations  of  the  Chinese  Digital  Seismographic  Network  (CDSN). 
Nguyen  and  Hsu  (1993)  inverted  surface  wave  group  velocity  dispersion  curves  from  the 
1^2  Yellow  Sea  earthquake  to  obtain  models  for  each  of  the  CDSN  stations.  Cipar  (1995) 
refined  their  focal  mechanism  using  a  grid  search  procedure  and  proposed  a  new  crustal 
structure  model  for  the  Yellow  Sea. 

The  1992  Yellow  Sea,  China,  Earthquake 

A  moderate  earthquake  occurred  in  the  Yellow  Sea  off  the  coast  of  eastern  China  on  3 
November  1992  at  17h  31m  23.7s.  Coordinates  are  35.3280N,  123.3120E;  focal  depth  is  10 
km,  mb  =  4.8  (Preliminary  Determination  of  Epicenters,  November,  1992,  Monthly  Listing). 
Figure  1  shows  the  earthquake  location  and  the  locations  of  the  CDSN  stations  used  in  this 
study  (Table  1).  Each  station  recorded  the  data  in  three  pass  bands:  long-period  (LH),  broad¬ 
band  (BH),  and  short-period  (SH)  (Peterson  and  TUgner,  1985).  Nguyen  (1994)  used  CDSN 
surface  wave  observations  to  determine  the  mechanism  of  the  1992  Yellow  Sea  earthquake  to 
be  nearly  strike-slip  with  a  minor  amount  of  dip-slip.  His  mechanism,  shown  in  Figure  1,  is 
consistent  with  first-motion  measurements  observed  on  CDSN  short-period  records.  The  fault 
plane  has  the  parameters:  strike  =  15°,  dip  =  80®,  and  rake  =  155®.  He  determined  the  source 
depth  and  seismic  moment  to  be  9  km  and  8.42x10^2  dyne-cm,  respectively. 

In  this  report,  I  use  the  CDSN  broadband  data  (the  BH  band)  to  place  an  additional 
constraint  on  the  upper  mantle  velocity.  The  data  set  is  less  than  ideal.  The  earthquake  was 
of  moderate  magnitude  located  at  sea  far  from  the  nearest  seismic  station,  hence  its  location  is 
imprecisely  known  and  teleseismic  observations  are  poor.  In  the  context  of  nuclear  treaty 

Table  1.  CDSN  Station  Information 


Station 

Location 

Latitude 

(ON) 

Lcxigitude 

(®E) 

Elevation 

(m) 

Delta 

(deg) 

Distance 

(km) 

Azimuth 

(deg) 

Back 

Azimuth 

(deg) 

BJI 

40.0403 

116.1750 

43.0 

7.357 

817.958 

311.860 

127.490 

Beiiing 

ENH 

30.2800 

109.4975 

487.0 

12.661 

1408.070 

250.471 

62.954 

Enshi 

HIA 

49.2667 

119.7417 

610.0 

14.175 

1575.558 

350.415 

167.994 

HaUar 

HKC 

22.3036 

114.1719 

0.0 

15.238 

1694.948 

214.049 

29.608 

Hong 

Kong 

KMI 

25.1500 

102.7500 

1952.0 

20.421 

2271.258 

245.854 

55.375 

Kunming 

LZH 

36.0867 

103.8444 

1560.0 

15.836 

1760.627 

278.427 

86.992 

Lanzhou 

MDJ 

44.6164 

129.5919 

250.0 

10.449 

1161.622 

25.507 

209.557 

Mudanjians 

SSE 

31.0956 

121.1867 

10.0 

4.581 

509.467 

203.475 

22.310 

Sheshan 

monitoring,  however,  it  is  precisely  these  weak  signals  which  will  be  of  prime  importance  for 
discrimination  of  suspect  events. 
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Figure  1.  Map  of  the  location  of  the  Yellow  Sea  earthquake  (star)  and  stations  (filled  triangles) 
used  in  this  study.  Open  triangles  denote  other  stations  in  the  area.  The  inset  shows  the  focal 
mechanism  reported  by  Nguyen  (1994);  shaded  quadrants  are  compressional. 
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Model  for  Crustal  Structure 

The  starting  velocity  model  for  the  crust  was  determined  by  Nguyen  and  Hsu  (1993) 
as  their  MDJ  model.  The  model  was  extended  into  the  upper  mantle  using  a  constant  velocity 
gradient.  Although  the  attenuation  was  set  rather  high,  Ae  values  are  broadly  consistent  with 
previous  measurements  for  Ae  continental  crust  (Cheng  and  Mitchell,  1981).  This  starting 
model  is  designated  A5  and  is  given  m  Table  2a. 

Cipar  (1995)  used  model  A5  and  Ae  long-period  (LH)  channels  at  CDSN  stations 
SSE,  MDJ,  BJl,  HIA,  LZH,  and  ENH  to  determme  Ae  focal  mechanism  by  a  gird  search 
Aversion  procedure.  The  best  fittmg  depA  is  15  km,  6  km  deeper  Aan  Ae  depA  determmed 
by  Nguyen  (1994).  The  inferred  mechanism  is  slightly  different  as  well  wiA  a  strike  =  19^, 
dip  =  90®,  and  rake  =  173®.  This  mechanism  is  broadly  consistent  wiA  Ae  results  reported 
by  Nguyen  (1994)  especially  given  Ae  fact  Aat  this  Aversion  did  not  use  any  constraAts  such 
as  P-wave  polarities.  It  should  be  noted,  however,  Aat  Ae  focal  mechamsm  reported  by 
Nguyen  (1994)  produced  raAer  large  misfits  at  all  source  depAs.  The  Aferred  source  depA, 
15  km,  is  well  withA  Ae  depAs  expected  for  a  contAental  earAquake. 

The  focal  mechamsm  grid  search  provides  a  source  model  Aat  produces  overall  better 
agreement  between  observed  and  synAetic  seismograms.  FurAer  improvement  can  be  gaAed 
by  modifyAg  Ae  crustal  structure.  Cipar  (1995)  tried  several  simple  modifications  of  Ae 
MDJ  paA  model  (model  A5).  The  most  successful  was  dubbed  A8  (Table  2b)  A  which  Ae 
lower  crustal  velocity  is  reduced  from  7.16  km/sec  for  A5  to  6.8  km/sec,  and  Ae  Q  reduced. 
The  Pn  velocity  of  AS  seems  to  be  quite  low  compared  to  oAer  continental  regions,  however, 
synAetics  computed  for  Wgher  Pn  velocity  values  produced  larger  misfits  Aan  A8.  Cipar 
(1995)  adopted  Ae  Nguyen  and  Hsu  (1993)  Pn  velocity  as  correct. 

Upper  Mantle  Velocity 

The  broadband  (BH  band)  data  for  Ae  Yellow  Sea  earAquake  are  shown  A  record 
section  format  A  Figure  2.  The  seismograms  have  been  Agh-pass  filtered  to  remove  long- 
period  noise  and  are  plotted  as  trace  normalized.  PromAent  phases  are  a  clear  Pn,  especially 
at  close  distances,  and  Lg.  A  straight  lAe  fit  to  Ae  first  arrivA  observations  at  SSE,  B JI, 

ENH,  and  HIA  (see  Figure  1),  over  Ae  distance  range  509  to  1575  km,  Adicates  Ae  average 
Pn  velocity  is  8.06  km/sec. 

Broadband  vertical  seismograms  for  Ae  two  closest  CDSN  stations  (SSE  and  BJI)  are 
shown  A  Ae  left  panel  of  Figure  3,  along  wiA  travel  time  curves  computed  for  model  A8 
(Table  2b).  We  are  lookAg  at  Ae  first  60  seconds  of  record,  Aat  is,  Ae  P  wave  group.  At 
these  ranges,  Ae  first  arrival  is  a  strong,  clear  Pn.  The  travel  time  curves  were  computed  for  a 
surface  source,  whereas  Ae  earAquake  occurred  at  10  to  15  km  depA.  NeverAeless,  Ae 
arrival  times  of  Ae  observed  Pn  phases  closely  match  Ae  Aeoretical  curves  in  time.  The  A8 
synAetics  (right  panel)  were  computed  usAg  Ae  reflectivity  meAod  (Fuchs  and  Muller, 

1971)  for  a  15-km  deep  source,  hence  Ae  time  discrepancy  between  travel  time  curves  and 
s)mAetics.  This  implies  Aat  Ae  travel  times  of  Ae  synAetics  are  early  compared  to  Ae  data. 
One  obvious  explanation  is  uncertaAty  A  location  and  origA  time.  Alternatively,  model  A8 
may  be  too  fast  relative  to  Ae  true  structure. 

The  synAetic  waveforms  provide  a  further  constraAt  on  Ae  strucAre.  One  of  Ae 
featores  of  model  A8  is  Wgh  velocity  gradients  A  Ae  upper  mantle  which  turn  energy  back  to 
Ae  surface,  leadAg  to  strong  first  arrivals  compared  to  later  phases.  The  data,  on  Ae  oAer 
hand,  show  strong  first  arrivals  followed  by  a  long  coda,  puncAated  by  several  Astinct 
phases.  Unlike  Ae  A8  synAetics,  however,  Ae  observed  Pn  phases  are  comparable  A 
ampliAde  to  later  phases. 

A  range  of  models  wiA  constant  gradients  A  Ae  upper  160  km  were  tested  by 
computing  synAetics  and  comparAg  Aese  to  Ae  observed  seismograms.  SynAetics  for 
model  A20  (Table  2c)  wiA  a  gradient  0.0013  km/sec/km  are  shown  A  Figure  4.  The  Pn 
ampliAde  at  BJI  is  comparable  qualitatively  to  Ae  observed  record.  The  synAetic  Pn  at  SSE, 
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Figure  2.  Broadband  CDSN  seismograms  of  the  1992  Yellow  Sea  earthquake  displayed 
as  trace  normalized  records.  The  seismograms  have  been  high-passed  filtered  using  a  4-pole 
Butterworth  filter  with  a  comer  frequency  of  0.4  Hz.  The  line  marked  8.06  km/sec  indicates  the 
best  fit  straight  line  to  measured  first  arrival  times. 
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on  the  other  hand,  is  too  small  compared  to  the  observed  Pn.  The  original  MDJ  model  (A5) 
has  a  Moho  velocity  of  7.76  km/sec  and  we  have  already  measured  the  average  Pn  velocity  to 
be  8.06  km/sec.  These  observations  suggest  a  strong  gradient  just  below  the  Moho,  which 
decreases  with  depth.  Note  that 


Table  2.  Crustal  Models 


Depth 

km 

P-velocity 

km/s 

Qp 

S-velocity 

km/s 

Qs 

Density 

g/cm^ 

0.0 

4.06 

(a) 

200.0000 

Model  A5 

2.3500 

100.0000 

2.3300 

1.0 

4.06 

200.0000 

2.3500 

100.0000 

2.3300 

1.0 

5.42 

1000.0000 

3.1300 

500.0000 

2.5800 

3.0 

5.42 

1000.0000 

3.1300 

500.0000 

2.5800 

3.0 

6.16 

1000.0000 

3.5600 

500.0000 

2.7500 

16.0 

6.16 

1000.0000 

3.5600 

500.0000 

2.7500 

16.0 

7.16 

1000.0000 

4.1400 

500.0000 

3.0300 

40.0 

7.16 

1000.0000 

4.1400 

500.0000 

3.0300 

40.0 

7.76 

1000.0000 

4.4800 

500.0000 

3.2300 

70.0 

8.10 

1000.0000 

4.6767 

500.0000 

3.3203 

100.0 

8.20 

1000.0000 

4.7344 

500.0000 

3.3582 

130.0 

8.40 

1000.0000 

4.8499 

500.0000 

3.4339 

160.0 

8.50 

1000.0000 

4.9076 

500.0000 

3.4718 

0.0 

4.0600 

(b)  Model  A8 

50.0000  2.3500 

25.0000 

2.3300 

1.0 

4.0600 

50.0000 

2.3500 

25.0000 

2.3300 

1.0 

5.4200 

100.0000 

3.1300 

50.0000 

2.5800 

3.0 

5.4200 

100.0000 

3.1300 

50.0000 

2.5800 

3.0 

6.1600 

525.0000 

3.5600 

263.0000 

2.7500 

16.0 

6.1600 

525.0000 

3.5600 

263.0000 

2.7500 

16.0 

6.8000 

525.0000 

3.9261 

263.0000 

2.8278 

40.0 

6.8000 

525.0000 

3.9261 

263.0000 

2.8278 

40.0 

7.7600 

525.0000 

4.4800 

263.0000 

3.2300 

70.0 

8.1000 

525.0000 

4.6767 

263.0000 

3.3203 

100.0 

8.2000 

525.0000 

4.7344 

263.0000 

3.3582 

130.0 

8.4000 

525.0000 

4.8499 

263.0000 

3.4339 

160.0 

8.5000 

525.0000 

4.9076 

263.0000 

3.4718 

40.0 

8.0000 

(c)  Model  A20 
(same  as  A8  above  40  km) 
525.0000  4.6200 

263.0000 

3.2300 

70.0 

8.0400 

525.0000 

4.6200 

263.0000 

3.3203 

100.0 

8.0800 

525.0000 

4.6300 

263.0000 

3.3582 

130.0 

8.1200 

525.0000 

4.6400 

263.0000 

3.4339 

160.0 

8.1600 

525.0000 

4.6415 

263.0000 

3.4718 
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Yellow  Sea  BH  Data  +  YSA-A8  Synthetics  &  Travel-Time  Curves 
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Figure  3.  Broadband  data  of  the  Yellow  Sea  earthquake  (left  panel)  and  synthetic  seismograms  computed 
for  model  A8.  Travel  time  curves  are  calculated  for  model  A8  at  zero  source  depth. 
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Figure  4.  Broadband  data  of  the  Yellow  Sea  earthquake  (left  panel)  and  synthetic  seismograms  computed 
for  model  A20.  Travel  time  curves  are  calculated  for  model  A20. 


170 


2(K«JL-B5 


the  travel  time  curves  for  model  A20  nearly  overlay  the  observed  Pn  first  arrivals. 
CONCLUSIONS  AND  RECOMMENDATIONS: 

The  average  upper  mantle  Pn  velocity  beneath  the  YeUow  Sea  and  eastern  China  is 
8.06  km/sec,  determined  using  data  from  the  1992  Yellow  Sea  earthquake,  and  is  consistent 
with  other  stable  platform  regions  (Meissner,  1986).  The  P  wave  velocity  gradient  is  0.0013 
km/sec/km,  although  it  may  be  higher  in  the  uppermost  mantle  just  below  the  Moho.  Array 
observations  of  a  suite  of  events  at  a  range  of  distances  and  azimuths  could  be  invaluable  in 
putting  constraints  on  crust  and  upper  mantle  velocities. 
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